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Abstract We present an application of the stereoscopic self-similar-expansion 
model (SSSEM) to Solar Terrestrial Relations Observatory (STEREO)/5'Mn- 
Earth Connection Coronal and Heliospheric Investigation (SECCHI) observa¬ 
tions of the 03 April 2010 CME and its associated shock. The aim is to verify 
whether CME-driven shock parameters can be inferred from the analysis of j- 
maps. Eor this purpose we use the SSSEM to derive the CME and the shock 
kinematics. Arrival times and speeds, inferred assuming either propagation at 
constant speed or with uniform deceleration, show good agreement with Ad¬ 
vanced Composition Explorer (ACE) measurements. The shock standoff distance 
[A], the density compression [^] and the Mach number [M] are calculated 
combining the results obtained for the CME and shock kinematics with models 
for the shock location. Their values are extrapolated to Li and compared to 
in-situ data. The in-situ standoff distance is obtained from ACE solar-wind 
measurements, and the Mach number and compression ratio are provided by the 
Harvard-Smithsonian Center for Astrophysics interplanetary shock database. 
They are ^ = 2.84 and M = 2.2. The best fit to observations is obtained 
when the SSSEM half width A = 40° and the CME and shock propagate with 
uniform deceleration. In this case we hnd A = 23R0, ^ = 2.61, and M = 2.93. 
The study shows that CME-driven shock parameters can be estimated from the 
analysis of time-elongation plots and can be used to predict their in-situ values. 
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1. Introduction 

Coronal mass ejections (CME) are large-scale eruptions of plasma and magnetic- 
field that originate from the Sun and propagate in the interplanetary medium. 
They arise above photospheric bipolar regions and travel outwards with speeds 
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between a few hundreds and more than 2500 km s (Yashiro et al. 


2004 


Cre- 


mades and Bothmer 2004). When the speed of the CME in the rest frame of 


the solar-wind exceeds the local Alfven speed, a magnetohydrodynamics (MHD) 
shock is driven ahead of the CME leading edge. Turbulence is generated in the 
post-shock region, the shock sheath, resulting in fluctuations of the magnetic- 
field and plasma parameters (Gosling and McComas 1987). Already in the 1970s, 
by comparing CME speeds at a few solar radii above the Sun to the Alfven speed 
at those heights, it was realized that CMEs can drive shocks in the corona (see, 

and references therein). With the help of 


Vourlidas and Ontiveros 2009 


e-g, _ 

today’s advanced coronagraphs, under suitable circumstances, these shocks can 
be detected in white-light images due to the density compression ahead of the 
CME (Vourlidas, 20061. Although indirect evidence of the presence of shocks can 
be identified, for example, in metric Type-II radio bursts and in observations of 
deflected ambient streamers, it was not until 2003 that the first direct detection 
of a CME-driven shock was announced. Vourlidas et al. (2003) employed an MHD 
model of the observed event, and they were able to show that the bright feature 
in the coronagraph images indeed corresponds to a shock in their simulation. 
Ontiveros and Vourlidas (2009) subsequently reported the analysis of LASCO 
fast events {v > 1500 km s“l), showing that in 86 % of the cases CME-driven 
shocks could be detected in the white-light images. Remote sensing observations 
of shocks can be exploited to determine shock parameters, such as compression 
ratio and Mach number (Maloney and Gallagher, 2011), and even further to 


estimate the coronal magnetic-field strength (Gopalswamy and Yashiro 2011 


Poomvises et al.. 2012). These estimates rely on early models for location of 


the Earth’s bow shock derived by Spreiter, Summers, and Alksne (19661 and 
Seiff (1962), and subsequently modified and extended to the low-Mach-number 
regime (Farris and Russell, 1994 Russell and Mulligan 2002). Gopalswamy and 
Yashiro (2011) and Poomvises et al. (2012) derived the quantities necessary for 
the calculation of the shock parameters via forward modeling using respectively 
a circular geometry and the graduated-cylindrical-shell (GCS) model by Th- 
ernisien, Howard, and Vourlidas (2006). 

In the present work, we employ inverse-modeling techniques to verify whether 
they allow the determination of the heliospheric-shock parameters and a predic¬ 
tion of their in-situ values. For this purpose we analyze white-light observations 
from the Sun Earth Connection Coronal and Heliospheric Investigation (SEC- 
CHI) suite onboard the NASA Solar Terrestrial Relations Observatory mission 
(STEREO: [Kaiser et al . 2008 Howard et al.. 


and its associated shock (Mostl et al. 2010 


2008) of the 03 April 2010 CME 


Liu et al. 2011). We separately 


determine the CME and shock kinematics and the time evolution of the standoff 
distance, the compression ratio, and the Mach number. Their values are then 
extrapolated to Li, where they can be compared to in-situ Advanced Composi¬ 
tion Explorer measurements (ACE: [Stone et oT 1998). 

The data processing and analysis techniques employed are presented in Section 
[2] Section [3| describes the models used in this work for the CME and shock 
kinematics and for the calculation of the shock parameters; their application to 
the 03 April 2010 CME is illustrated in Section [^ In Section [^ we discuss the 
applicability of the SSSEM for the derivation of the shock parameters. Section 
[^provides a summary and conclusions. 
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2. Data and Image Processing 


The NASA/STEREO mission, launched on 25 October 2006, consists of the 
twin spacecraft, STEREO-A (Ahead) and -B (Behind) orbiting the Sun at 1 
AU. STEREO-A precedes the Earth in its orbit, while STEREO-B follows it. 
The separation angle between each spacecraft and Earth increases by about 22.5° 
per year. A package of five telescopes, the SECCHI suite, is mounted onboard 
both STEREO-A and -B, and provides observations of the solar corona from the 
solar disk to distances beyond 1 AU (Howard et al, 2008). 

The EUVI instrument images the Sun’s atmosphere in four different wavelengths, 
from the upper chromosphere to the corona; CORl and COR2 are two nested 
coronagraphs with fields of view covering respectively the ranges from 1.5 to 
4 R 0 , and from 2.5 to 15R0. HIl and HI2 are heliospheric imagers observing 
the inner heliosphere up to distances larger than 1 AU. HIl has a field of 20°, 
centered at 14° from the Sun along the Sun-Earth line; HI2 has a field of view 
of 70°, centered at 54° along the Sun-Earth line. 

Combined observations from the instruments of the SECCHI suite allow us to 
follow the evolution of solar-wind features from their origin near the solar surface 
up to beyond 1 AU (Davis et al. 2009). Multi-point observations, moreover, pro¬ 
vide an unprecedented insight on the three-dimensional structure and evolution 
of CMEs, and allow the determination of their kinematics free from projection 
effects. 

In order to use STEREO data for scientific analysis, processing of the images 
is required. Elat-held and vignetting correction, CCD bias subtraction, and 
correction for optical distortion are some of the calibration procedures applied to 
COR2 observations (for a complete description of COR images processing visit 
hesperia.gsfc.nasa.gov/ssw/stereo/secchi/doc/cor_prep.html). For HIl and HI2 
observations, background stray-light removal is necessary to reveal the presence 


of CMEs (Davies et al. 2009 Eyles et al. 2009). 


The visibility of faint features in white-light observations can be enhanced by 
image differencing, a technique in which each image is obtained by subtracting 


two subsequent snapshots (Sheeley et al. 1999). What is shown in difference 


images is the change in brightness between two consecutive snapshots: white 
features lead black ones in the direction of motion, the former corresponding to 
plasma which, in the later image, has propagated further away from the Sun. 
Time-elongation plots (also known as j-maps) are built by stacking slices of 
difference images cut at a fixed position angle, and they can be used to deter¬ 


mine the kinematics of CMEs (Davies et al. 2009); structures propagating in 


the imagers’ field of view appear as curved tracks in time-elongation plots (see 
Eigurej^. In the interpretation of Sheeley et al. (1999) the shape of such curves 
is due to the apparent acceleration/deceleration seen by an observer measuring 
the angular distance (i.e. the elongation [e]) of a feature that propagates away 
from the Sun at constant speed. 
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Figure 1. Time-elongation plot for the time interval between 02 April 2010 and 07 April 2010 
obtained from STEREO-A COR2, HIl and HI2 difference images. The track corresponding to 
the leading edge (LE) of the CME analyzed in this work is the one appearing in the COR2 
field of view on 03 April 2010 at about 09:24UT. In the elongation range observed by HIl a 
faint curve appears, preceding the CME. This corresponds to the density enhancement due to 
compression downstream of the shock. The tracks following the CME leading edge correspond 
to regions of density enhancement in the CME itself (also see Figure!^. Manual selection 
of points along the track corresponding to the CME and the shocked plasma allows for the 
independent determination of their kinematics. 


3. Methods 


In order to determine the CME and shock kinematics we build j-maps from cuts 
of STEREO-A and -B difference images along the Ecliptic plane (see Figure 
[^. This allows us to determine the temporal evolution of the CME and shock 
height, speed, and direction of propagation in the Ecliptic plane. The results 
can therefore be directly compared to ACE measurements. Manual selection of 
points along the curved tracks visible in the j-maps yields the evolution of the 
features’ elongation angles as a function of time. The time-elongation profile 
thus obtained can be transformed to distances from the Sun provided a set 
of assumptions on the CME geometry. In this article we use the stereoscopic 


self-similar expansion model (SSSEM) developed by Davies et al. (2012 2013), 
which represents the CME front as a circular front expanding at a fixed angular 
width [A] (Figure [^. Via geometrical arguments it is possible to show that at 
each instant of time, for each spacecraft [i], the distance from the Sun of the 
propagating feature is related to the measured elongation via 
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R = di 


sin(ei)(l + sin(A)) 
sin(ei + 4>i) + sin(A) ’ 


i = A,B, 


( 1 ) 


with di being the Sun-spacecraft distance, [e^] the elongation angle obtained 
from the j-maps, and (j)i being the direction of propagation of the white-light 
feature with respect to the Sun-spacecraft line, while A is an input parameter 
of the model. With the additional constraint 


7AB — a + IpB 


( 2 ) 


where jab is the (known) separation angle between the two spacecraft, si¬ 
multaneous STEREO-A and -B observations allow us to directly solve Equations 
Q for R, tpA and (/>b and derive the kinematic properties of solar-wind features. 
The direction of propagation with respect to the Sun-Earth line is simply 
derived as the difference between the Sun-spacecraft angle [Ji] and the direction 
of propagation with respect to the Sun-spacecraft line determined from the 
SSSEM 

The speed is obtained, using a quadratic Lagrangian interpolation, as the nu¬ 


merical time derivative of the radial profile (Davies et al, 20131. It should be 


noted that Equation is only valid when STEREO-A, STEREO-B the CME, 
Sun, and Earth lie in the same plane. This condition is satished in our study 
since we follow the propagation of the CME and the shock in the Ecliptic plane. 
The model has two limiting cases: for A = 0° it coincides with the fixed-((> 


approximation (Sheeley et al. 1999 Rouillard et al. 2008), which represents the 
CME as a point source. As in Equation ([^, (/) is the direction of propagation 
with respect to the Sun-observer line. The case A = 90°, instead, corresponds 
to the harmonic mean geometry, describing the CME as a sphere anchored at 


the Sun (Lugaz et al., 2010). 

It is then possible, from time series of white-light observations, to determine the 
temporal evolution of outward moving features in the spacecraft held of view. 
If the CME and the shock are visible and distinguishable in white-light images 
and time-elongation plots, application of the SSSEM allows to separately derive 
the CME and shock kinematics. We applied the SSSEM for different values of 
A to the time-elongation prohles manually derived from the j-maps. We studied 
the interval between A = 0° and A = 90° varying A in steps of 10°. 

The errors in the determined height [ar] and direction of propagation are 

obtained using the approach of Liu et al. (2010b). These are calculated via the 
following relations: 


= 



( 3 ) 





d(j)A 

deA 


<^A 


+ 



( 4 ) 
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Following the approach of Liu et at. (2010bI, we assign an uncertainty of ten 
pixels to the measured elongation angles, corresponding to 0.04° for points in 
COR2, 0.2° for points in HIl, and 0.7° for points in HI2. The derived errors are 
a function of the direction of propagation [4>i\, the height, the half-width A, the 
positions of the spacecraft, and the instruments’ resolution. The errors in the 
speed and in the shock parameters are obtained by propagating the errors in R 
and (f)i. 

In order to infer the in-situ arrival time and speed we fit the height-time profile to 
a linear function. The method also yields the eruption time. The fit is restricted 
to points in the HIl field of view, where both the CME and the shock velocity can 
be approximated as constant. A discussion on the accuracy of constant-velocity 
fits to the observed elongation profiles has been given by Mbstl et al. (2014). 
The comparison between arrival time and speed prediction based on remote¬ 
sensing observations and in-situ measurements shows that, independently of 
the fitting method (FP, HM, SSEM), the predicted arrival speeds are higher 
than the observed ones, and the corresponding arrival times are earlier. This is 
supported by Lugaz and Kintner (2013), who argue that solar-wind drag, which 
has the effect of decelerating fast CMEs towards the solar-wind speed, has a 
substantial impact on the accuracy of arrival-time predictions. Eor this reason 
we also perform a quadratic fit to the height profile in the HIl field of view, to 
verify whether the inclusion of deceleration improves the accuracy of the in-situ 
extrapolations. Arrival times and speeds are corrected accounting for the CME 
and shock geometry and direction of propagation following Mostl and Davies 
(2013). When the CME or the shock hit ACE, at a distance d = 213R0 from 
the Sun, the position of their apex is determined from the angle of propagation 
[5 b] and the half width [A] according to: 


Rsssem = d 


1 -|- sin(A) 


cos{5e) + \/sin(A)2 — sin(5B)^ 


(5) 


The importance of the correction increases with increasing A and 6e- We 
remind the reader that A S [0°,90°] and 6e G [—180°, 180°], with Se = 0° 
corresponding to propagation along the Sun-Earth line. 

The upstream Mach number and the compression ratio can be derived from 
the standoff distance, i.e. the distance between the shock front and the CME 
leading edge, by employing models initially developed for the Earth’s bow shock. 


Seiff (1962) and Spreiter, Summers, and Alksne (1966), based on hydrodynamic 


theory and experimental data, showed that for the Earth’s bow shock, in the 
high-Mach-number regime, the standoff distance primarily depends on the ratio 
between upstream and downstream solar-wind densities [^] according to 


^ = 1 . 1 -. 

D pd 


( 6 ) 


A is the standoff distance and D is the position of the nose of the magne¬ 
tosphere measured from Earth, i.e the size of the dayside magnetosphere in a 


frame of reference with its origin at Earth (Farris and Russell, 1994). 
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Using hydrodynamic relations (Landau and Lifshitz, 1959) it can be shown that 
the density compression ratio [^] can be related to the upstream Mach number 
via 


Pu 

Pd 


(7 - 1) M2 + 2 


( 7 ) 


(7+1) M2 ’ 

with 7 being the adiabatic index. We adopted the value 7 = 5/3, as in Farris 


and Russell (1994). 


Farris and Russell (1994) modified Spreiter’s relation, replacing the size [D] in 
Equation ®> with the radius of curvature 


A = 0.81^, 

Pd 


( 8 ) 


arguing that the shock location would be more sensitive to the local shape of 
the magnetosphere, represented by the radius of curvature, rather than to its 
size alone. The change in the proportionality parameter between Equations ®> 
and (1^ results from the fact that Farris and Russell (|1994|) estimated for the 


magnetosphere = 1.35. Furthermore, they introduced a —1 term in the 
denominator on the right hand side of Equation ([^ , 


A (7-1)M2 + 2 

R, (7 + 1)(M2-1)’ 


( 9 ) 


extending the applicability of Equations ([^ and ([^ to the low-Mach-number 
regime: the presence of the —1 term, in fact, allows the shock to move to infinity 
as the Mach number approaches unity. 

Although a full MHD approach would be required to completely describe the 
CME-shock system, as well as the interaction between the solar-wind and the 
Earth’s magnetosphere, it has been shown that the above approach yields good 
results in the high-Mach-number regime, for M > 5 (Spreiter, Summers, and 


Alksne, 1966). It is also a good approximation in the low-Mach-number regime, 


provided that M is substituted for the fast magnetosonic Mach number. When 
the Alfven speed substantially exceeds the sound speed, moreover, the magne¬ 
tosonic Mach number is not too different from the Alfvenic Mach number [Ma] 
(Fairfield et al., 2001), and M can be replaced by Ma- Equation © has been 
used to derive CME and shock properties based on in-situ (Russell and Mul- 


ligan, 

2002 

) and remote-sensing observations 

'Gopalswamy and Yashiro 

2011 

Poomvises et al. 

2012|Maloney and Gallagher 

2011 

). The analysis of white-light 


images, in particular, allows the determination of shock parameters, and even of 
the magnetic-field strength, at distances from the Sun where no direct plasma 
measurements are available. In this work we derive the standoff distance as the 
difference between the shock and the CME height determined via the SSSEM. 
The SSSEM, moreover, yields the CME radius of curvature as a function of the 
CME half-width [A], according to 
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Figure 2. Schematic representation of the stereoscopic self-similar-expansion model for 
A = 40° (left panel) and its limiting cases: the fixed-0 (top right panel) and the harmonic mean 
(bottom right panel). According to the SSSEM the CME is approximated to a sphere expanding 
with a constant half width [A]. The positions of the STEREO spacecraft in the pictures match 
the location of STEREO-A and -B at the time of observations. The elongations ca a^nd ee are 
obtained from time—elongation plots for STEREO-A and -B observations. The Sun-spacecraft 
distances [dp^] and [^b] as well as the separation between the spacecraft, 7 ab = are 

known at all times. Given the half width [A] as an input parameter, it is possible to determine 
the CME height above the Sun [R] and its direction of propagation with respect to STEREO-A 
and -B, [0 a and 0 b] as a function of time. The fixed-0 and harmonic mean approximations 
are the limiting cases of the SSSEM, for A = 0° and A = 90° respectively. 


Rc 


= R 


sin(A) 

1 + sin(A) ’ 


( 10 ) 


The temporal evolution of the compression ratio and the Mach number can 
therefore be derived using Equations Q and 0 via an analysis of j-maps.By 
assuming propagation at constant velocity or constant acceleration it is possible 
to extrapolate these quantities to the position of the ACE spacecraft, and com¬ 
pare them with in-situ measurements. We remind the reader that, in order to 
determine the standoff distance at Li, we have to take the difference between 
the shock and CME height corrected according to Equation 0 - 


4. Results 


The event analyzed in this work erupted from NOAA Active Region 11059 (S25° 
W03°) on 03 April 2010. The CME eruption was associated with a B7.4 flare 
detected by GOES at 09:04 UT, and resulted in a geomagnetic storm with Kp 
8 - on 05 April 2010 ( Mostl et aZ.[|2010{|Liu et al. 20111. The separation angles 
between the STEREO-A and -B spacecraft and Earth at the time of observations 
are 67° and 71° respectively (Figure [^. 

In-situ measurements of the CME and its associated shock are provided by the 
ACE spacecraft: the shock passage, marked by the sudden increase in velocity 
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Figure 3. a) STEREO-A (left panel) and -B (right panel) observations of the 03 April 2010 
CME. The compressed region ahead of the CME leading edge is visible throughout the entire 
CME propagation, and particularly clear in HIl images. The Milky Way in the HI2B field of 
view compromises the analysis of stereoscopic observations at the elongations imaged by the 
instrument. 


and density, is detected on 05 April 2010 at about 08:00UT. It is followed after 
about four hours by the CME. Its signatures can be identified, for instance, 
in the decreasing velocity profile, in the variation of the (j) component of the 
magnetic-field, associated to a smooth rotation of the field itself, and in the 
lower temperature compared to the solar-wind (Figure]^. 

Figurej^shows difference images obtained from STEREO-A and -B COR2, HIl, 
and HI2 observations. HI2B images are subject to high levels of background 
noise due to the transit of the Milky Way through the instrument’s field of view. 
Stereoscopic observations of the transient are therefore only available up to 24 
degrees in elongation, corresponding to the outer field of view of the HIl imager. 
The CME morphology seems to be preserved during its outward propagation. 
Although the images in Figure are orthogonal to the Ecliptic plane, where we 


SOLA: Volpes_Bothiner2015.tex; 11 September 2015; 0:27; p. 9 







10 


L. Volpes, V.Bothmer 



Figure 4. ACE plasma and magnetic-field measurements for the time interval 04 April 2010 
to 07 April 2010. The shock passage is marked by the sudden increase in density and velocity 
at 08:00, UT on 05 April 2010; it was followed by the CME, which crossed the spacecraft at 
12:00 and lasted for about 25 hours. The narrow gray area indicates the shock sheath and the 
wider pink area corresponds to the CME. 

The vertical lines represent the CME and shock-arrival times estimated via a second-degree 
fit to the stereoscopic self-similar expansion model. The light and dark dashed green lines 
mark the CME arrival for A = 0° and A = 90°, respectively. The light black dot-dashed line 
corresponds to the shock-arrival time for A = 0°, and the dark black dot—dashed line to the 
shock arrival for A = 90°. 


determine the CME and shock kinematics, the absence of significant distortion 
is a good indication that self-similar-expansion is a reasonable assumption for 
the event under study. This is because what we observe in white-light is the 
integrated election density along the line of sight. 

A region of enhanced density ahead of the CME, corresponding to compressed 
plasma downstream of the shock, can be detected in difference images of the 
event. It is especially bright in HIl observations, and it remains visible even in 
the near-Earth environment imaged by HI2 A. 

Figurej^shows the j-map for the event analyzed in this work. It was obtained by 
juxtaposing cuts of STEREO-A images along the Ecliptic plane for observations 
obtained between 02 April 2010 and 07 April 2010. The elongation range in the 
y axis covers the fields of view of COR2, which extends up to 4°, HIl and HI2, 
respectively between about 4° and 24°, and 19° and 89°. 
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Figure 5. Temporal evolution of the height, direction of propagation, and velocity for the 03 
April 2010 CME (left panel) and the shock (right panel), derived from an application of the 
stereoscopic self-similar-expansion model to STEREO-A and -B COR2 and HIl observations. 
The orange dots indicate solutions of the fixed-^ approximation (A = 0°); the black dots 
represent the solution to the harmonic mean approximation (A = 90°). The CME and the 
shock appear to propagate radially in the Ecliptic plane from the CME’s source region, located 
at S25°W03° on the Sun. 


The bright, extended curve appearing on 03 April 2010 at about 09:24 UT 
follows the evolution of the CME leading edge. A second, fainter curve preceding 
the CME can be identified in elongation range observed by the HI instrument. 
It corresponds to the propagation of the faint feature ahead of the CME visible 
in HIl images (see Eigurej^, which we associate with compressed plasma down¬ 
stream of the shock. Since the two tracks are distinguishable from one another 
in the j-maps obtained from both STEREO-A and -B observations (the latter 
is not shown), the CME and shock kinematics could be separately determined. 
This was done by applying the SSSEM to the time-elongation profile derived in 
the range corresponding to the COR2 and HIl fields of view. 

Eigure illustrates the CME and shock height, speed, and direction of propa¬ 
gation with respect to Earth for the cases A = 0° and A = 90°. Solutions for 
the intermediate values of A lie between the two sets of points in each plot, and 
are therefore not shown. The average direction of propagation of the CME with 
respect to the Sun-Earth line, Se = 5° ± 2° for A = 0° and = 6° ± 2° for 
A = 90°, is compatible with radial evolution (in the Ecliptic plane) from its 
source region, located at W03° latitude on the solar disk. These results do not 
differ much from the findings of Mostl et al. (20101, who determined 6e = 9° for 
A = 0° and = — 5° for A = 90°, and = 0° ± 2° from geometrical forward 
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Table 1. Upper panel: CME arrival and eruption time, and shock-arrival time 
and speed determined via a constant-speed fit to the self-similar-expansion model 
(SSSEM) results for the HIl fields of view. ACE in-situ measurements of the shock 
and CME arrival are reported for comparison, as well as the time of flaring activity. 
In both the cases A = 0° (FP approximation) and A = 90° (HM approximation) the 
derived eruption times are very close to the observed flare time, with a difference of 
about 30 minutes for the FP results and 1 hour for the HM results. Predicted arrival 
times are roughly eight hours earlier than the time of CME and shock passage at the 
location of ACE, and consistently the speeds are higher. Lower panel: Same quantities 
as above, obtained via a constant-negative-acceleration fit to the SSSEM results. The 
eruption times are again in excellent agreement with observations, and a significant 
improvement of the arrival-time and speed prediction is found, with accuracy as high 
as two hours. The inferred speeds are also closer to the ones determined through the 
linear fit. 



ACE-I-EUVI 

1st degree fit A = 0° 

1st degree fit A = 90° 

CME eruption time 
Shock arrival time 
Shock arrival speed 
CME arrival time 

03 Apr. 2010 09:04 
05 Apr. 2010 07:56 
750 km s~^ 

03 Apr. 2010 12:24 

03 Apr. 2010 08:24 

05 Apr. 2010 00:26 

1025 km s~^ 

05 Apr. 2010 05:24 

03 Apr. 2010 08:32 

05 Apr. 2010 00:19 

1033 km s~^ 

05 Apr. 2010 06:05 



ACE-I-EUVI 

2nd degree fit A = 0° 

2nd degree fit A = 90° 

CME eruption time 
Shock arrival time 
Shock arrival speed 
CME arrival time 

03 Apr. 2010 09:04 
05 Apr. 2010 07:56 
750 km s~^ 

03 Apr. 2010 12:24 

03 Apr. 2010 08:27 

05 Apr. 2010 03:38 

867 km s~^ 

05 Apr. 2010 08:44 

03 Apr. 2010 08:12 

05 Apr. 2010 04:19 

851 km s~^ 

05 Apr. 2010 10:48 


modeling. 

Eruption time and arrival time and speed at the location of ACE are inferred 
by applying a first- and a second-degree fit to the CME and shock heights in 
the range imaged by HIl only. Table shows a comparison between fits and 
observations. Eruption times show an excellent agreement with the time of 
the GOES fiare detection and the time of the CME first appearance in the 
CORl field of view, at 09:15 UT for both the linear and the quadratic fit. The 
overestimation of the propagation speed for both the CME and the shock in the 
constant-velocity fit leads to estimated arrival times earlier than the observed 
ACE crossing. Including deceleration significantly improves the arrival-time and 
speed predictions, particularly for the A = 90° case. Based on the CME and shock 
heights obtained from the SSSEM results, the standofi" distance, the compression 
ratio, and the Mach number can be computed according to Equations (§, 
and (10), substituting 7 = 5/3. Their time profiles and the evolution of the 
standoff distance and the radius of curvature are shown in Figure for A = 0° 
and A = 90°. 

The case A = 0° is excluded, as the corresponding zero curvature leads to 
a singularity in Equation ( [^. It should be pointed out that the compression 
ratio [y] plotted in Figure |6l) is the inverse of the ratio on the left-hand side 
of Equations (|^ and ([^. Both the density compression and the Mach number 
appear to be roughly constant during the CME evolution, in agreement with 
the results found by Maloney and Gallagher (2011 1 . The increase of the standoff 
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Date and Time (UT) 


L 



3-ADr-2010 3-Apr-2010 3-Apr-2010 3-Apr-2010 

06:40 12:13 17:46 23:20 

Date and Time (UT) 



3-Apr-2010 3-/^r-2010 3-Apr-2010 3-Apr-2010 

f2:29 15:49 19:09 22:29 


Date and Time (UT) 


Date and Time (UT) 


Figure 6. Upper panel: standoff distance and CME radius of curvature as a function 
of time obtained via the application of the stereoscopic self-similar-expansion model 
(SSSEM) to STEREO-A and -B COR2 and HIl observations, for A = 10° (orange dot), 
and A = 90° (light black diamonds). Lower panel: temporal evolution of the compression 
ratio and Mach number derived according to Equations ||^ and using the values of 
the standoff distance and the radius of curvature obtained from the SSSEM. The red dots 
represent SSSEM solutions for A = 10° and the black diamonds for A = 90°. The case A = 0° 
is not presented, as the zero radius of curvature would lead to a singularity in the Equation ©■ 


distance is also in good agreement with observations, which show the shock 
progressively moving away from its driver. 

An extrapolation of the standoff distance, the compression ratio, and the 
Mach number to the location of ACE was performed assuming again propagation 
at constant speed as well as at constant acceleration. The correction in Equation 
(§ was applied to the CME and shock heights to determine the standoff distance 
at Li. For the A = 10° case we find an unphysical compression ^ < 1; this can 
be caused by the overestimated value of the standoff distance, or the small CME 
curvature associated to this value of the half width. Validation is necessary for 
the extrapolated values of compression ratio and the Mach number at 1 AU. 
The in-situ standoff distance is calculated directly from the CME and shock- 
arrival time and speed as measured by ACE. The extrapolated values of the 
density compression ratio and the Mach number are instead compared to the 
results published on the interplanetary (IP) shock database of the Harvard- 
Smithsonian Center for Astrophysics (www.cfa.harvard.edu/shocks/). The shock 
parameters are determined applying different analysis methods and principles 
such as, e.g, the magnetic and velocity coplanarity theorem and the solution of 
the Rankine-Hugoniot equations, to the ACE and Wind plasma and magnetic- 
field measurements. We use the average of the results obtained with the different 
methods. 

Table shows the values of the standoff distance, Mach number, and density 
compression extrapolated to 1 AU via linear and quadratic fit to the SSSEM 
results, as well as those calculated from in-situ data. The in-situ standoff dis¬ 
tance is Aace = 19R0. A comparison with the linear and quadratic fits to 
the SSSEM results, shows that for all the values of A the extrapolated standoff 
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Table 2. Upper panel: extrapolation to 1 AU of the standoff distance, the 
Mach number, and the compression ratio. The results have been obtained 
via first-degree fit to the CME and shock height-time profile determined 
from the self-similar-expansion model (SSSEM) applied to STEREO-A 
and -B HI images, for A = 10° and A = 90°. The value of the standoff 
distance determined from in-situ ACE measurements and those of the 
Mach number and compression ratio available on the Harvard-Smithso- 
nian Center for Astrophysics interplanetary shock database are reported 
for comparison. Lower panel: same quantities as above, determined via 
a second-degree fit to the CME and shock height. The standoff distance 
derived via the SSSEM is higher than the observed one. The values of 
the compression ratio and the Mach number computed from in-situ mea¬ 
surements are within the range defined by our results. For A = 10°, the 
unphysical result ^ < 1 can be related to the high standoff distance or 
to the small associated CME curvature. 



ACE 

1st degree fit A = 

10° 

1st degree fit A ; 

- 90° 

Standoff distance 

19 R0 

30 Rq 


26 Rq 


Mach number 

2.2 

1.38 


3.26 


Compression ratio 

2.84 

0.75 


2.82 




ACE 

2nd degree fit A = 

^ 10° 

2nd degree fit A 

0 

o 

05 

II 

Standoff distance 

19 Rq 

29 Rq 


26 Rq 


Mach number 

2.2 

1.41 


3.46 


Compression ratio 

2.84 

0.78 


2.91 



distance, ranging from 2 IR 0 to SOR©, overestimates the observed one. This 
is easily recognizable in Figure where the measured and the extrapolated 
magnetosheath thicknesses can be directly compared. 

The in-situ value of the Mach number, Mace = 2.2 falls within the ranges deter¬ 
mined via the linear and quadratic ht, 1.38 < M < 3.50 and 1.40 < M < 3.54, 
respectively. The same is true for the compression ratio, whose in-situ value 
— = 2.84 is contained in the intervals derived from first and second-degree 

fits, 0.75 < < 2.97, 0.78 < < 2.91. The value of A that best approximates 

the in-situ Mach number and compression is A = 40° for both the linear and 
quadratic fit to the SSSEM height. For the linear fit we find A = 27R0, ^ = 2.19 
and M = 2.42, while for the quadratic ht we hnd A = 23R0, ^ = 2.61 and 
M = 2.93. As for the arrival time we hnd a signihcant improvement in the in-situ 
extrapolations when deceleration is included in our hts. 


5. On the Errors Associated to OfF-Ecliptic Propagation 

In order to derive the Mach number and the compression ratio using the relations 
in Equations ([^ and ([^, the standoff distance has to be measured at the CME 
nose. Therefore attention should be paid to the direction of propagation of the 
CME when employing the SSSEM for the determination of shock parameters. 
If the direction of propagation of the CME with respect to the Ecliptic plane 
[a] is not small the determination of shock parameters using the SSSEM might 
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Figure 7. Graduated-cylindrical-shell fitting for the 03 April 2010 CME for the COR2 and HIl 
fields of view. The fit performed in COR2 was extended to HIl assuming self-similar-expansion. 
The CME propagates at constant longitude and latitude. The latter was used as input for the 
assessment of the errors associated to the calculation of the shock parameters from the standoff 
distance computed via the SSSEM. 


lead to incorrect results. The standoff distance would be measured at the wrong 
location and the SSSEM radius of curvature might not well represent the one at 
the nose of the CME. Furthermore, in order to compare the SSSEM results to in- 
situ data, we computed the standoff distance from the SSSEM heights corrected 
according to Equation ([^ . In this way we determined the shock parameters from 
the standoff distance calculated along the Sun-Earth line. It is now necessary to 
estimate how much these values deviate from those computed at the CME nose. 
From STEREO observations (see Figure]^ it is evident that the CME propagates 
southwards both in COR2 and in HIl. To estimate its direction of propagation 
in 3D we model the CME and the shock with the graduated-cylindrical-shell 
(CCS) model by Thernisien, Howard, and Vourlidas (20061. The model, which 
describes the CME as a croissant-shaped flux rope expanding self-similarly, has 
previously been employed to determine the shock parameters (Poomvises et al. 
(2012)). We start from results obtained in COR2, and extend them to the HIl 
fields of view. We expand the CCS flux rope keeping the aspect ratio and the 
half angle constant, i.e. assuming self-similar propagation. We find that the 
CME propa_gates southwards at an angle of a = —23° with the Ecliptic plane 
(see Figure!^, and at an angle of about = 5° with respect to the Sun- 

Earth line. This is in agreement with our SSSEM results, which yield Se ~ 1° 
for the shock and 6e ^ 5° for the CME. No deflection is observed during the 
CME propagation, i.e. the CME longitude and latitude [a] are the same 
in the COR2 and HI fields of view. We use the value of a determined from the 
CCS model and that of 5^; obtained from the SSSEM to estimate the difference 
between the standoff distance and the CME curvature computed at the nose and 
the ones derived from the SSSEM. We assume a to be the same for the CME 
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and the shock. Finally, we determine the Mach number and compression ratio 
at the CME nose and compare them to the results in Table 
To estimate the value of the standoff distance at the CME nose we assume that 
the CME is a sphere of radius i?, as this is the straightforward extension of 
the SSSEM to three-dimensional space. We assume that the intersection of the 
sphere with the Ecliptic plane is the SSSEM circle (see Figure [Sa]). We now 
consider the plane formed by the CME height vector (ie. the vector originating 
at the Sun and ending at Rsssem) and perpendicular to the ecliptic. The plane 
cuts the sphere at the poles, so that the projection of the sphere on this is a 


circle of radius R (see Figure 8b). The SSSEM circle in this plane degenerates 


to a segment of length 2r, where r is the radius of the SSSEM circle, i.e. 


r = Rsssem 


sin(A) 

1 -I- sin(A) 


( 11 ) 


The radius of the sphere, which corresponds to the “true” CME curvature, is 
given in terms of r according to 


R = 


sin(/3) ’ 


( 12 ) 


(see Figure [8b|). 

The angle /3 can be determined considering that BC = AB tan /3, with AB = r, 
and that BC = OB tana, with OB = Rsssem — r. We get 


tan/3 = 


[Rsssem — r) tana 
Finally, the CME apex is at a distance H = OC + R given by 


(13) 


H = 


[Rsssem - r) 


sin /3 


(14) 


At a fixed latitude, the error in neglecting the off-Ecliptic propagation will 
decrease with increasing A, i.e. with increasing size of the CME. This is because 


the horizontal axis in Figure 8b will slice the sphere closer to its Equator for 
larger CMEs. At fixed A, on the other hand, the error will increase with increasing 
a, as the slice will be cut closer to the poles, i.e. further away from the apex. 

It is possible then to derive the standoff distance and the CME curvature at 
the CME nose, compute the shock parameters, and compare them to our in-situ 
extrapolations. The results are summarized in Table for the case A = 40°, 
which corresponds to our best fit to the data. For each parameter X in the 
Table [i.e. standoff distance, curvature, compression ratio, and Mach number) 
we report the value of , A^orr being the parameter calculated at the CME 

apex. 

The radius of curvature is the quantity that is affected the most by neglecting 
the off-Ecliptic propagation of the CME and the shock. The apex value of the 
standoff distance is at most 5 % bigger than the one we determined with the 
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(a) (b) 

Figure 8. Left panel: CME propagation along the Ecliptic plane, as in Figure]^ Right panel: 
CME propagation perpendicular to the Ecliptic plane, with the horizontal axis aligned to the 
CME height vector. The plane intersects the CME, assumed to be a sphere, at its poles. 

Table 3. Relative difference in the values of the shock parameters 
calculated along the Sun-Earth line and at the CME apex for the 
best fit to in-situ observations obtained via application of the 
SSSEM. 


Relative error for: 

1st degree fit A = 40° 

2nd degree fit A = 40° 

Standoff distance 

5 % 

4 % 

Curvature 

20 % 

20 % 

Mach number 

13 % 

23 % 

Compression ratio 

13 % 

15 % 


SSSEM. The difference between the SSSEM values of the Mach number and 
compression ratio computed at the CME nose and via the SSSEM are between 
23 % and 13 %, with the second-degree fit being affected the most. This indicates 
that the results we obtained by computing the standoff distance via the SSSEM 
do not deviate dramatically from those obtained calculating the standoff distance 
and the curvature of the CME at its nose. 


6. Summary and Conclusions 


It has been long believed that fast CMEs can drive shocks in the corona, and 
that their signatures could be observed in coronagraph images due to the density 
compression across the shock. Thanks to the increasing instrument sensitivity 
it is now possible to detect CME-driven shocks in a vast number of events 
(Ontiveros and Vourlidas, 2009). Stereoscopic observations, moreover, provide 
an unprecedented insight into the 3D properties of coronal mass ejections, and 
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allow the determination of CME kinematics without the need to assume constant 


velocity and strictly radial propagation (Liu et al. , 2010a 2010b). In this work we 
present the analysis of the fast CME that erupted on 03 April 2010, associated 
with an interplanetary shock detected by ACE on 05 April 2010. The results of 
this work show that the shock was already present in the corona below ISR©, and 
was visible in white-light images throughout the whole evolution towards Earth. 
Exploiting the fact that the CME and the shock are visible and distinguishable 
in time-elongation plots, we derived the CME and shock kinematics indepen¬ 
dently, employing the SSSEM for different values of the angular width. To our 
knowledge, this is the first time that the analysis of j-maps has been applied 
separately to a CME and its driven shock. The results of our study show that 
the CME evolved radially in the Ecliptic plane from the eruption site located at 
S25° W03° on the solar disk. Linear and second-degree fits to the height profile 
were performed, showing that the latter yields better estimates for arrival times 
and speeds. The standoff distance, the compression ratio, and the Mach number 
evolution were derived; their values were extrapolated to the location of ACE 
by exploiting a linear and a quadratic fit to the CME and shock heights. In 
order to validate these results, the standoff distance was computed based on 
the in-situ CME and shock arrival times and speeds. The values of the Mach 
number and the compression ratio were compared to results published on the 
CfA IP shock database, computed from ACE plasma and magnetic-field data, 
according to which M = 2.2 and ^ = 2.84. We found the case A = 40° to best 
match observations. If constant velocity is assumed, A = 27Rq, ^ =2.19 and 
M = 2.42. When deceleration is assumed to take place in the CME propagation, 
we find A = 23R0, ^ = 2.61 and M = 2.93. 

Although the CME apex was at an angle a = —23° with respect to the Ecliptic 
plane, an assessment of the error associated to the use of the SSSEM to deter¬ 
mine the shock parameters showed that their values at the CME nose do not 
substantially differ from those obtained via the SSSEM. 

The methods presented in this work allow the determination of CME-driven 
shock properties from remote-sensing observations, where in-situ measurements 
are not available. Magnetic-field strength estimates are also possible based on 
such results (Gopalswamy and Yashiro, 20111, which overcome the difficulty of 
directly measuring the coronal magnetic-field. 

The study shows that it is possible to use inverse-modeling techniques to deter¬ 
mine the time evolution of CME-driven shock parameters and to predict their 
in-situ values. The NASA Solar Probe Plus mission to be launched in 2018 will 
explore the near Sun environment at distances up to 9 . 86 R 0 , i.e. at 0.046 AU, 
and the ESA Solar Orbiter mission, also to be launched in 2018, will orbit the 
Sun at a distance of 0.28AU. The white-light observations from the Wide-Field 
Imager on Solar PRobe Plus (WISPR: Vourlidas et al. 20151 and Solar Orbiter 
Heliospheric Imager (SoloHI: Howard et al. 2013) cameras coupled with in-situ 
measurements will provide unprecedented knowledge about the coronal plasma 
and magnetic-field properties at locations unexplored so far. Application of the 
models employed in this work to CME images acquired by WISPR and SoloHI 
will give insight on the spatial variation of CME-driven shock properties. At the 
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same time, a deeper knowledge of the plasma and magnetic-field characteristics 
at the same locations will help to further validate the model results. 

Acknowledgments The authors acknowledge funding by Deutsche Forschungsgemeinschaft 
(DEG) under grant SFB 963/1, project A2. They are further grateful for support in the analysis 
of STEREO SECCHI data by Jackie Davies and Richard Harrison from Rutherford Appleton 
Laboratory, Didcot, Chilton, UK, for the information on shock databases provided by Luciano 
Rodriguez from the Royal Observatory of Belgium, Bruxelles, and for the GCS modeling in 
COR2 by Eckhard Bosman from the University of Gottingen. Volker Bothmer acknowledges 
support of the CGAUSS (Coronagraphic German And US Solar Probe Plus Survey) project 
for WISPR by the German Space Agency DLR under grant 50 OL 1201. 


References 

Cremades, H., Bothmer, V.: 2004, On the three-dimensional configuration of coronal mass 
ejections. Astron. Astrophys. 422, 307. DOI ADS 

Davies, J.A., Harrison, R.A., Rouillard, A.P., Sheeley, N.R., Perry, C.H., Bewsher, D., 
Davis, C.J., Eyles, C.J., Crothers, S.R., Brown, D.S.: 2009, A synoptic view of solar 
transient evolution in the inner heliosphere using the Heliospheric Imagers on STEREO. 
Geophys. Res. Lett. 36, 2102. DOI ADS 

Davies, J.A., Harrison, R.A., Perry, C.H., Mdstl, C., Lugaz, N., Rollett, T., Davis, C.J., 
Crothers, S.R., Temmer, M., Eyles, C.J., Savani, N.P.: 2012, A Self-similar Expansion Model 
for Use in Solar Wind Transient Propagation Studies. Astrophys. J. 750, 23. DOI ADS 

Davies, J.A., Perry, C.H., Trines, R.M.G.M., Harrison, R.A., Lugaz, N., Mdstl, C., Liu, Y.D., 
Steed, K.: 2013, Establishing a Stereoscopic Technique for Determining the Kinematic Prop¬ 
erties of Solar Wind Transients based on a Generalized Self-similarly Expanding Circular 
Geometry. Astrophys. J. 777, 167. DOI ADS 

Davis, C.J., Davies, J.A., Lockwood, M., Rouillard, A.P., Eyles, C.J., Harrison, R.A.: 2009, 
Stereoscopic imaging of an Earth-impacting solar coronal mass ejection: A major milestone 
for the STEREO mission. Geophys. Res. Lett. 36, 8102. DOI ADS 

Eyles, C.J., Harrison, R.A., Davis, C.J., Waltham, N.R., Shaughnessy, B.M., Mapson-Menard, 
H.C.A., Bewsher, D., Crothers, S.R., Davies, J.A., Simnett, G.M., Howard, R.A., Moses, 
J.D., Newmark, J.S., Socker, D.G., Halain, J.-P., Defise, J.-M., Mazy, E., Rochus, P.: 2009, 
The Heliospheric Imagers Onboard the STEREO Mission. Solar Phys. 254, 387. DOI ADS 

Fairfield, D.H., Cairns, I.H., Desch, M.D., Szabo, A., Lazarus, A.J., Aellig, M.R.: 2001, The 
location of low Mach number bow shocks at Earth. J. Geophys. Res. 106, 25361. DOI 
ADS 

Farris, M.H., Russell, C.T.: 1994, Determining the standoff distance of the bow shock: Mach 
number dependence and use of models. J. Geophys. Res. 99, 17681. DOI ADS 

Gopalswamy, N., Yashiro, S.: 2011, The Strength and Radial Profile of the Coronal Magnetic 
Field from the Standoff Distance of a Coronal Mass Ejection-driven Shock. Astrophys. J. 
Lett. 736, L17. DOI ADS 

Gosling, J.T., McComas, D.J.: 1987, Field line draping about fast coronal mass ejecta - A 
source of strong out-of-the-ecliptic interplanetary magnetic fields. Geophys. Res. Lett. 14, 
355. DOI ADS 

Howard, R.A., Moses, J.D., Vourlidas, A., Newmark, J.S., Socker, D.G., Plunkett, S.P., Ko- 
rendyke, C.M., Cook, J.W., Hurley, A., Davila, J.M., Thompson, W.T., St Cyr, O.C., 
Mentzell, E., Mehalick, K., Lemen, J.R., Wuelser, J.P., Duncan, D.W., Tarbell, T.D., 
Wolfson, C.J., Moore, A., Harrison, R.A., Waltham, N.R., Lang, J., Davis, C.J., Eyles, 
C.J., Mapson-Menard, H., Simnett, G.M., Halain, J.P., Defise, J.M., Mazy, E., Rochus, P., 
Mercier, R., Ravet, M.F., Delmotte, F., Auchere, F., Delaboudiniere, J.P., Bothmer, V., 
Deutsch, W., Wang, D., Rich, N., Cooper, S., Stephens, V., Maahs, G., Baugh, R., Mc- 
Mullin, D., Carter, T.: 2008, Sun Earth Connection Coronal and Heliospheric Investigation 
(SECCHI). Space Sci. Rev. 136, 67. DOI ADS 

Howard, R.A., Vourlidas, A., Korendyke, C.M., Plunkett, S.P., Carter, M.T., Wang, D., Rich, 
N., McMullin, D.R., Lynch, S., Thurn, A., Clifford, G., Socker, D.G., Thernisien, A.F., 
Chua, D., Linton, M.G., Keller, D., Janesick, J.R., Tower, J., Grygon, M., Hagood, R., 
Bast, W., Liewer, P.C., DeJong, E.M., Velli, M.M.C., Mikic, Z., Bothmer, V., Rochus, P., 


SOLA: Volpes_Bothiner2015.tex; 11 September 2015; 0:27; p. 19 



20 


L. Volpes, V.Bothmer 


Halain, J.-P., Lamy, P.L.: 2013, The solar and heliospheric imager (SoloHI) instrument for 
the solar orbiter mission. In: Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series 8862, 0. DOI ADS 

Kaiser, M.L., Kucera, T.A., Davila, J.M., St. Cyr, O.C., Guhathakurta, M., Christian, E.: 

2008, The STEREO Mission: An Introduction. Space Sci. Rev. 136, 5. DOI ADS 
Landau, L.D., Lifshitz, E.M.: 1959, Fluid mechanics, Pergamon Press, Oxford. ADS 
Liu, Y., Davies, J.A., Luhmann, J.G., Vourlidas, A., Bale, S.D., Lin, R.P.: 2010a, Geometric 
Triangulation of Imaging Observations to Track Coronal Mass Ejections Continuously Out 
to 1 AU. Astrophys. J. Lett. 710, L82. DOI ADS 
Liu, Y., Thernisien, A., Luhmann, J.G., Vourlidas, A., Davies, J.A., Lin, R.P., Bale, S.D.: 
2010b, Reconstructing Coronal Mass Ejections with Coordinated Imaging and in Situ Ob¬ 
servations: Global Structure, Kinematics, and Implications for Space Weather Forecasting. 
Astrophys. J. 722, 1762. DOI ADS 

Liu, Y., Luhmann, J.G., Bale, S.D., Lin, R.P.: 2011, Solar Source and Heliospheric Conse¬ 
quences of the 2010 April 3 Coronal Mass Ejection: A Comprehensive View. Astrophys. J. 
734, 84. DOI ADS 

Lugaz, N., Kintner, P.: 2013, Effect of Solar Wind Drag on the Determination of the Properties 
of Coronal Mass Ejections from Heliospheric Images. Solar Phys. 285, 281. DOI ADS 
Lugaz, N., Hernandez-Charpak, J.N., Roussev, LI., Davis, C.J., Vourlidas, A., Davies, J.A.: 
2010, Determining the Azimuthal Properties of Coronal Mass Ejections from Multi- 
Spacecraft Remote-Sensing Observations with STEREO SECCHI. Astrophys. J. 715, 493. 
DOI ADS 

Maloney, S.A., Gallagher, P.T.: 2011, STEREO Direct Imaging of a Coronal Mass Ejection- 
driven Shock to 0.5 AU. Astrophys. J. Lett. 736, L5. DOI ADS 
Mostl, C., Davies, J.A.: 2013, Speeds and Arrival Times of Solar Transients Approximated by 
Self-similar Expanding Circular Fronts. Solar Phys. 285, 411. DOI ADS 
Mostl, C., Temmer, M., Rollett, T., Farrugia, C.J., Liu, Y., Veronig, A.M., Leitner, M., Galvin, 
A.B., Biernat, H.K.: 2010, STEREO and Wind observations of a fast ICME flank triggering 
a prolonged geomagnetic storm on 5-7 April 2010. Geophys. Res. LeR. 37, 24103. DOI ADS 
Mostl, C., Amla, K., Hall, J.R., Liewer, P.C., De Jong, E.M., Colaninno, R.C., Veronig, A.M., 
Rollett, T., Temmer, M., Peinhart, V., Davies, J.A., Lugaz, N., Liu, Y.D., Farrugia, C.J., 
Luhmann, J.G., Vrsnak, B., Harrison, R.A., Galvin, A.B.: 2014, Connecting Speeds, Direc¬ 
tions and Arrival Times of 22 Coronal Mass Ejections from the Sun to 1 AU. Astrophys. J. 
787, 119. DOI ADS 

Ontiveros, V., Vourlidas, A.: 2009, Quantitative Measurements of Coronal Mass Ejection- 
Driven Shocks from LASCO Observations. Astrophys. J. 693, 267. DOI ADS 
Poomvises, W., Gopalswamy, N., Yashiro, S., Kwon, R.-Y., Olmedo, O.: 2012, Determination 
of the Heliospheric Radial Magnetic Field from the Standoff Distance of a CME-driven 
Shock Observed by the STEREO Spacecraft. Astrophys. J. 758, 118. DOI ADS 
Rouillard, A.P., Davies, J.A., Forsyth, R.J., Rees, A., Davis, C.J., Harrison, R.A., Lockwood, 
M., Bewsher, D., Crothers, S.R., Eyles, C.J., Hapgood, M., Perry, C.H.: 2008, First imaging 
of corotating interaction regions using the STEREO spacecraft. Geophys. Res. Lett. 35, 
10110. DOI ADS 

Russell, C.T., Mulligan, T.: 2002, On the magnetosheath thicknesses of interplanetary coronal 
mass ejections. Planet. Space Sci. 50, 527. DOI ADS 
Seiff, A.: 1962, Gas Dynamics in Space Exploration. In: The NASA University Conference on 
the Science and Technology of Space Exploration, NASA SP 24, NASA, Washington DC, 
19. 

Sheeley, N.R., Walters, J.H., Wang, Y.-M., Howard, R.A.: 1999, Continuous tracking of coronal 
outflows: Two kinds of coronal mass ejections. J. Geophys. Res. 104, 24739. DOI ADS 
Spreiter, J.R., Summers, A.L., Alksne, A.Y.: 1966, Hydromagnetic flow around the magneto¬ 
sphere. Planet. Space Sci. 14, 223. DOI ADS 

Stone, E.C., Frandsen, A.M., Mewaldt, R.A., Christian, E.R., Margolies, D., Ormes, J.F., 
Snow, F.: 1998, The Advanced Composition Explorer. Space Sci. Rev. 86, 1. DOI ADS 
Thernisien, A.F.R., Howard, R.A., Vourlidas, A.: 2006, Modeling of Flux Rope Coronal Mass 
Ejections. Astrophys. J. 652, 763. DOI ADS 

Vourlidas, A.: 2006, Detections of CME-Driven Shocks with LASCO. In: Lacoste, H., Ouwe- 
hand, L. (eds.) SOHO-17. 10 Years of SOHO and Beyond, ESA SP SP-617, ESA, 
Noordwijk, 23. ADS 


SOLA: Volpes_Bothiner2015.tex; 11 September 2015; 0:27; p. 20 



An Application of the SSSEM to the Determination of CME-Driven Shock Parameters 21 


Vourlidas, A., Ontiveros, V.: 2009, A Review of Coronagraphic Observations of Shocks Driven 
by Coronal Mass Ejections. In: Ao, X., Burrows, R., Zank, G. (eds.) Shock Waves in space 
and astrophysical environments, AIP CS 1183, AIP, ???, 139. DOI 
Vourlidas, A., Wu, S.T., Wang, A.H., Subramanian, P., Howard, R.A.: 2003, Direct Detection of 
a Coronal Mass Ejection-Associated Shock in Large Angle and Spectrometric Coronagraph 
Experiment White-Light Images. Astrophys. J. 598, 1392. 

Vourlidas, A., Howard, R.A., Plunkett, S.P., Korendyke, C.M., Thernisien, A.F.R., Wang, D., 
Rich, N., Carter, M.T., Chua, D.H., Socker, D.G., Linton, M.G., Morrill, J.S., Lynch, S., 
Thurn, A., Van Duyne, P., Hagood, R., Clifford, G., Grey, P.J., Velli, M., Liewer, P.C., Hall, 
J.R., DeJong, E.M., Mikic, Z., Rochus, P., Mazy, E., Bothmer, V., Rodmann, J.: 2015, The 
Wide-Field Imager for Solar Probe Plus (WISPR). Space Sci. Rev.. DOI ADS 
Yashiro, S., Gopalswamy, N., Michalek, G., St. Cyr, O.C., Plunkett, S.P., Rich, N.B., Howard, 
R.A.: 2004, A catalog of white light coronal mass ejections observed by the SOHO spacecraft. 
J. Geophys. Res. (Space Phys.) 109, 7105. 


SOLA: Volpes_Bothiner2015.tex; 11 September 2015; 0:27; p. 21 



